Astronomy & Astrophysics manuscript no. 7050-sissa 


© ESO 2008 


February 1, 2008 





o 
o 

>> 



6 



> 

O 

in 
o 
i> 
o 



X 



A new comprehensive set of elemental abundances in DLAs 

III. Star formation histories ^ 

M. Dessauges-Zavadsky\ F. Calura^, J. X. Prochaska^, S. D'Odorico'^, and F. Matteucci^'^ 

' Observatoire de Geneve, Universite de Geneve, 51 Ch. des Maillettes, 1290 Sauvemy, Switzerland 

- INAF - Osservatorio Astronomico di Trieste, Via G. B. Tiepolo 11, 34131 Trieste, Italy 

^ UCO/Lick Observatory, University of California at Santa Cruz, Santa Cruz, CA 95064, USA 

* European Southern Observatory, Karl-Schwarzschildstr. 2, 85748 Garching bei MUnchen, Germany 

^ Dipartimento di Astronomia, Universita di Trieste, via G.B. Tiepolo 11, 34131 Trieste, Italy 

Received; accepted 

ABSTRACT 

We obtained comprehensive sets of elemental abundances for eleven damped Lyor systems (DLAs) at zdla = 1.7-2.5. For nine of 
them, we accurately constrained their intrinsic abundance patterns accounting for dust depletion and ionization effects. In Paper I of 
this series, we showed for three DLA galaxies that we can derive their star formation histories and ages from a detailed comparison 
of their intrinsic abundance patterns with chemical evolution models. We determine in this paper the star formation properties of six 
additional DLA galaxies. The derived results confirm that no single star formation history explains the diverse sets of abundance 
patterns in DLAs. We demonstrate that the various star formation histories reproducing the DLA abundance patterns are typical of 
local irregular, dwarf starburst and quiescent spiral galaxies. Independent of the star formation history, the DLAs have a common 
characteristic of being weak star forming galaxies; models with high star formation efficiencies are ruled out. The distribution of the 
DLA star formation histories shows a trend of finding more galaxies with a star formation history typical of dwarf irregulars with a 
bursting star formation toward high redshifts, z > 2. Only two DLA galaxies (each at z < 2) in our sample of nine objects have a star 
formation history typical of spiral galaxies. Since DLAs sample the broad distribution of galaxies at high redshift, this trend indicates 
that young and less evolved proto-galactic structures with low masses and low star formation rates are more common toward higher 
redshifts. This is further supported by the star formation rate and age distributions. Indeed, all the derived DLA star formation rates 
per unit area are moderate or low, with values between -3.2 < logSFR < -1.1 Mq yr"' kpc"-. The DLA abundance patterns also 
require a large spread in ages ranging from 20 Myr up to 3 Gyr. Enhanced a over iron-peak ratios are associated with young objects 
having undergone a recent burst of star formation, while solar a over iron-peak ratios are associated with old objects undergoing an 
inefficient continuous star formation. The oldest DLA in our sample is observed at zdla = 1-864 with an age estimated to more than 
3 Gyr; it nicely indicates that galaxies were already forming at Zf > 10. But, most of the DLAs show ages much younger than that 
of the Universe at the epoch of observation. Young galaxies thus seem to populate the high redshift Universe at z > 2, suggesting 
relatively low redshifts of formation (z ~ 3) for most high-redshift galaxies. The large dispersion in star formation history and age 
indicates that the DLAs are drawn from a diverse population of galaxies. The DLA star formation properties are compared with those 
of other high-redshift galaxies identified in deep imaging surveys with the aim of obtaining a global picture of high-redshift objects. 
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1. Introduction The connection between the population of DLA absorbers 

and galaxy populations identified at high redshifts in deep imag- 
Quasar (QSO) absorption line systems provide an alterna- ing studies remains, however, a key topic in extragalactic astron- 
tive means to study the high-redshift galaxy properties other omy today. While an evolutionary link between the emission- 
than the traditional stellar luminosity selected galaxy samples selected galaxies at z = 2-3 and the bulge c omponent of 
dSteidel et al.ll2003l and reference therein). Indeed, they are ex- the galaxy population today is strongly suggeste d dSteidel et alj 



pected to select galaxies based on gas cross-section, rather than 1 19961; lAdelberger et al.ll2004i ISteidel et aLll2005h . an analogous 

luminosity, and hence give a complementary view of the galaxy connection for DLAs is still a matter of debate , because of 

population to that obtained from emission-selected galaxies. The the difficult detection of DLA galaxy counterparts (' M0ller et"al] 

damped Lya systems (DLAs) with H I column densities higher ,2004i) . Indeed, the complete picture relating the gas cross- 

than 2x10^'' cm"^ are particularly interesting, since they contain section, ionization and chemical abundances of the interstellar 

the majority of the neutral gas at high redshifts and may provide media and gaseous halos surrounding galaxies to the emission 

the fuel reservoirs for subsequent star for mation in galaxies (e.g. properties of the associated galaxy remains incomplete. 

IWolfe et aT]|1986l;|Prochaska et al.ll2005h . DLA galaxy counterparts are difficult to identify, because 

they typically lie at angular separations of the order of one arc- 
sec from the quasar line of sight, and are therefore swamped 



5e;jdc#pn;jfre9t(e.s?i'to:miroslava.dessauges@obs.unige.ch by the light from the quasar. At high redshift (z > 2), de- 

* Based on UVES observations made with the European Southern spite large observational efforts, emission, in general from Lya, 

Observatory VLT/Kueyen telescope, Paranal, Chile, collected during from seven DLA counterparts was detected and spectroscopi- 

the programme ID No. 70.B-0258(A). cally confirmed, but only four are not associated DLAs with 
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the Q SO, zdla ^ zqso (see a summary by IWeatherlev et alj 
12005 ). At lower redshift (z < 1), the DLA galaxy counterparts 
are more readily observed. Deep-imaging and follow-up spec- 
troscopic surveys of the fields of DLA absorbers revealed host 
galaxies with a wide range of morphological types between spi- 
rals, irregulars, dwarfs, and low surface brightness gal axies (e.g. 
iLe Brun et al]fT997l: lRaoltani200l Ichen et aLllIOOBl) . 

Studies of elemental abundances in high-redshift DLAs 
have, on the other hand, achieved large statistical samples over 
the last few years. This progress is the result of 8-10 m class 
telescopes and the high resolution spectroscopy. The HIRES 
spectrograph on the Keck I telescope at Mauna Kea, Hawaii, 
was the first in operation and allowed the routine detection of 
8 elements in DLAs: Al, Si, S, Cr, M n, Fe, Ni, and Zn (e.g. 
iLu et al.lll996l:lProchaska & Wolfdl 19991) . With the arrival of the 
UVES spectrograph on the Very Large Telescope (VLT) ESO 
telescope at Cerro Paranal, Chile, with a better efficiency in the 
blue (A < 4000 A) and in the red {A > 7000 A), * 7 addi- 
tional element abundances can now be accurately measured in 
DLAs: C, N, O, Mg, P, Ar, and Ti, plus the deuterium and the 
molecular hydrogen (e.g. iLevshakov et al.ll2002l: iPetitiean et alj 
120001) . Moreover, for some specific DLAs, the so-called metal- 
strong DLAs, the access to 10 new elements, several beyond 
the Fe-peak, is be coming very promising (Prochaska et a lj2003l: 
iHerbert-Fort et al]|2006h . 

The main motivation common to all these observational ef- 
forts is to understand the physical processes at play in the forma- 
tion and evolution of galaxies. To achieve this goal, the chem- 
ical information can be used as one of the possible means at 
our disposal to link the properties of DLA high-redshift galax- 
ies with both those of present-day galaxies and other galax- 
ies detected at high redshifts. Indeed, with the measure of so 
many diff'erent ele mental abundances, DLAs are r eal physi- 
cal laboratories. In |Dessauges^Z^adi^etal] ([2003, hereafter 
Paper I), we showed for three DLAs with comprehensive sets 
of elemental abundances (22 elements detected) that we can 
derive their star formation histories (SFHs) and ages from a 
detailed comparison of their abundance patterns with chemi- 
cal evolution models. Previous studies also investigated the na- 
ture of DLA galaxies by means of chemical evolution models 
(Jimenez, Bowen & Matteucci 1999: Hou et al. <2001l:ICora et alj 
[2003; Calura et al . 2003). Our approach diff'ers in the fact that 
we treat the DLAs studied as individual systems, and not as an 
ensemble following a common evolutionary sequence. 

We propose to constrain the nature of six new high redshift 
DLA galaxies via their chemical abundance analysis following 
the prescriptions defined in Paper I in order to derive a global 
picture of high-redshift galaxies. Together with the three DLAs 
previously studied, these nine DLA systems at zdla = 1.7 - 2.5 
were selected only on the basis of: (i) a bright background 
quasar, (ii) the existence of HIRES/Keck spectra, and (iii) mini- 
mal Lyo' forest contamination which maximizes access to chem- 
ical elements. The latter criterion was achieved by demanding 
(zqso - Zdla) < 0.6 while requiring that the velocity separation 
Av(QSO - DLA) > 3000 km s ' such that the two ai-e not physi- 
cally associated. They form a sample which covers the full range 
of H I column densities, metallicities, and dust depletion levels as 
observed in DLAs. The DLAs studied in this paper thus are well 
representative of the average high-redshift DLA galaxy popula- 
tion. 

The layout of the paper is as follows. In Sect. |2] we describe 
our DLA sample and the determination of their intrinsic abun- 
dances free from dust depletion and ionization effects. In Sect. [3] 



we present the chemical evolution models we will use to con- 
strain the DLA star formation histories and ages. In Sect. |4] we 
derive for the six DLA galaxies in our sample the star forma- 
tion histories and ages. Finally, in Sect. |5] we compare and try 
to connect the determined DLA properties, i.e. ages, star for- 
mation rates, and star formation histories, with the properties of 
emission-selected galaxies at the same redshifts. To conclude. 
Sect. |6]provides a summary of the main constraints obtained on 
the DLA galaxies. The adopted cosmology throughout the paper 
is Hq = 65 km s"' Mpc"' (h = 0.65), = 0.3, and Qa = 0.7. 



2. Our DLA sample and their abundances 

The sample of DLA galaxies that we will study in this paper 
is unique. Indeed, by combining our high quality UVES/VLT 
spectra with existing HIRES/Keck spectra, we obtained com- 
prehensive sets of elemental abundances for eleven DLAs. We 
detected 54 metal-line transitions which allowed us to mea- 
sure the column densities of 30 ions from 22 elements: B, C, 
N, O, Mg, Al, Si, P S, CI, Ar, Ti, Cr, Mn, Fe, Co, Ni, Cu, 
Zn, Ge, As, and Kr. This contrasts with the majority of DLAs 
for which only a ha ndful of elements (Si, Fe, Cr, Zn, Ni) is 
usually detected (e.g.lLuetai] |I996l: iProchaska & Wolfelll999l: 
iProchaska et al.ll200 lV In these other cases, one struggles to in- 
terpret the DLA abundance patterns, especially to distinguish 
between d ifferential depletion and in trinsic nucleosynthetic pat- 
terns (e.g. IProchaska & Wolfell2002l) . 

The abundance analysis of our sample of DLAs is 
presented in Paper I for four DLAs toward QBOlOO+13, 
Q B1331-I-17, QB2231-0 02, and QB2343+125, and in 
De ssauges-Zavadsky et al.l (|2006, hereafter Paper II) for seven 
DLAs toward QB0450-1310, QB0841-1-129, PKS 1157-1-014, 
QB1210+175, QB2230-i-02, and QB2348-14443. A complete 
interpretation of the abundances with the determination of the 
star formation histories can also be found in Paper I for the 
first set of DLAfl In this paper we will focus on the second 
set of DLAs. In Table [T] we summarize their ionic column 
densities and absolute abundances. The most important relative 
abundances are also reported. They have been computed by con- 
sidering only the column densities of Voigt profile components 
detected in both elements of the abundance ratio. 

The interpretation of the observed elemental abundance pat- 
terns in DLAs is not straightforward. The principal difficulty is to 
disentangle the nucleosynthesis contributions from dust deple- 
tion and ionization effects. Indeed, since we are measuring gas- 
phase elemental abundances, in presence of dust the observed 
abundances may not represent the intrinsic chemical composi- 
tion of the galaxy if part of the elements is removed from the gas 
to the solid ph ase, as it happens in the int erstellar medium (ISM) 
of our Galaxy dSavage & SembachI 19961) . Similarly, a fraction of 
the gas may be ionized by external or internal ionizing sources. 
Therefore, a comprehensive set of ionic column densities and 
elemental abundances is of great advantage in order to evaluate 
the dust depletion and ionization effects, and hence derive the 
intrinsic abundances. 



' The naming of quasars have slightly been modified relative to 
Papers I and II in order to be in accordance with the lAU rules. 

- The DLA toward Q B2343+ 125 was excluded from the star forma- 
tion history study, because of its large ionization fraction. 



Table 1. Summary of absolute and relative abundances 





Quasar 


QB()841 


+ 129 


QB0841 


+ 129 


PKS 1157+014 


QB12I0+175 


QB2230+02 


QB2348 


-1444 




Sola 


2.375 


2.476 




1.944 


1.892 


1.864 


2.279 


logWCHi) 


20.99(0.08) 


20.78(0.08) 


21.60(0.10) 


20.63(0.08) 


20.83(0.05) 


20.59(0.08) 


Ionic column densities and absolute abundances 


log ^(B n) 


|B/Hl„b, 










< 11.93 


< -0.46 














log N(N i) 


[N/HJ„b, 


14.60(0.01) 


-2.29(0.10) 


13.94(0.08) 


-2.76(0.13) 


> 15.29 


> -2.22 


14.71(0.09) 


-1.82(0.13) 


15.02(0.06) 


-1.68(0.11) 


13.35(0.05) 


-3.16(0.11) 


log A^(0 i) 


[0/H]„te 


> 16.05 


> -1.77 


16.15(0.10) 


-1.46(0.14) 














> 15.04 


> -2.38 


log N(Mg n) 


[Mg/H]„i,s 


15.13(0.10) 


-1.37(0.13) 


14.99(0.15) 


-1.35(0.17) 


15.99(0.04) 


-1.18(0.11) 


< 15.51 


< -0.53 


< 15.61 


< -0.80 


< 14.80 


< -1.37 




[Mg/H],-„(EO()) 




-1.37(0.13) 








-1.17(0.12) 




< -0.53 




< -0.77 




< -1.37 




[Mg/H]„,-(E1I) 




-1.37(0.13) 








-1.17(0.12) 




< -0.53 




< -0.77 




< -1.37 


log N(A\ u) 


|AI/H]„bs 


> 13.72 


> -1.76 


> 13.37 


> -1.90 


> 15.01 


> -1.08 


> 14.94 


> -0.18 






< 12.71 


< -2.37 


log ^(Si n) 


lSi/Hl„i„ 


15.21(0.04) 


-1.34(0.09) 


14.99(0.03) 


-1.35(0.09) 


15.97(0.01) 


-1.19(0.10) 


15.33(0.03) 


-0.86(0.09) 


15.70(0.03) 


-0.69(0.06) 


14.18(0.06) 


-1.97(0.10) 




[Si/H]^or(EO()) 




-1.34(0.09) 








-1.19(0.10) 




-0.86(0.09) 




-0.68(0.07) 




-1.96(0.10) 




[Si/H]„„(E11) 




-1.34(0.09) 








-1.18(0.11) 




-0.86(0.09) 




-0.68(0.07) 




-1.96(0.10) 


log N(P ii) 


[P/H]„bs 


12.82(0.06) 


-1.73(0.12) 


12.56(0.09) 


-1.78(0.13) 


13.86(0.09) 


-1.30(0.15) 


< 13.95 


< -0.24 


13.69(0.08) 


-0.70(0.11) 






log N(S ii) 


[S/H]„bs 


14.69(0.05) 


-1.48(0.11) 


14.48(0.10) 


-1.50(0.15) 






14.96(0.02) 


-0.85(0.10) 


15.29(0.06) 


-0.66(0.11) 


13.75(0.06) 


-2.03(0.12) 


log N(Cl i) 


[Cl/H]„bs 










< 12.81 


< -2.07 














log N(Ai i) 


[Ar/H]„bs 


13.53(0.08) 


-1.86(0.13) 


13.11(0.08) 


-2.07(0.13) 


















logWlTiu) 


[Ti/H]„bs 










12.82(0.13) 


-1.71(0.16) 


12.34(0.08) 


-1.06(0.13) 


12.68(0.09) 


-1.09(0.10) 






logW(Cru) 


[Cr/H]„bs 


13.07(0.02) 


-1.59(0.08) 


12.89(0.06) 


-1.58(0.1 1) 


13.76(0.02) 


-1.52(0.10) 


13.28(0.03) 


-1.02(0.09) 


13.44(0.04) 


-0.97(0.08) 


12.30(0.09) 


-1.97(0.12) 


logWCMnu) 


[Mn/H]„bs 


12.50(0.02) 


-2.02(0.08) 


12.33(0.15) 


-1.98(0.17) 


13.25(0.02) 


-1.87(0.10) 


12.72(0.01) 


-1.42(0.09) 


13.03(0.03) 


-1.33(0.06) 


11.66(0.15) 


-2.46(0.17) 


logiV(Fen) 


[Fe/H]„b. 


14.76(0.01) 


-1.73(0.08) 


14.50(0.03) 


-1.78(0.09) 


15.46(0.02) 


-1.64(0.10) 


15.01(0.03) 


-1.12(0.09) 


15.24(0.03) 


-1.09(0.06) 


13.84(0.05) 


-2.25(0.09) 




[Fe/H]c„(EOO) 




-1.59(0.10) 








-1.36(0.18) 




-0.98(0.16) 




-0.76(0.14) 




-2.13(0.19) 




[Fe/H],„(E11) 




-1.59(0.10) 








-1.36(0.18) 




-0.98(0.16) 




-0.75(0.14) 




-2.13(0.19) 


log N(Co II) 


[Co/H]„b. 










13.10(0.13) 


-1.40(0.16) 














logWlNiii) 


lNi/H]„b, 


13.53(0.05) 


-1.69(0.09) 


13.19(0.09) 


-1.84(0.13) 


14.21(0.02) 


-1.66(0.10) 


13.67(0.06) 


-1.19(0.11) 


14.09(0.04) 


-0.92(0.08) 


< 12.30 


< -2.30 




|Ni/H]„,(EOO) 




-1.56(0.11) 








-1.37(0.18) 




-1.03(0.19) 




-0.67(0.13) 




< -2.17 




|Ni/H]„,(ElI) 




-1.55(0.11) 








-1.38(0.18) 




-1.05(0.18) 




-0.59(0.16) 




< -2.18 


logWlCuii) 


|Cu/H]„bs 










< 12.31 


< -1.56 














log N(Zn II) 


lZn/H]„b, 


12.10(0.02) 


-1.49(0.09) 


11.69(0.10) 


-1.74(0.14) 


12.99(0.05) 


-1.27(0.12) 


12.40(0.05) 


-0.88(0.10) 


12.72(0.05) 


-0.67(0.09) 


< 11.28 


< -1.97 


logA'(Geii) 


[Ge/H]„b, 










11.99(0.14) 


-1.24(0.18) 














logWCAs II) 


[As/H]„b. 










< 11.78 


< -0.19 


< 11.97 


< +0.97 










logWCKri) 


[Kr/H]obs 










< 13.01 


< +0.10 


< 12.48 


< +0.54 










logWCCii*) 




12.96(0.08) 




< 13.08 




15.13(0.09) 








13.95(0.02) 




< 12.25 





Relative abundances 



[Zn/Fe]„b. 


+0.24(0.02) 


+0.04(0.11) 


+0.37(0.07) 


+0.24(0.07) 


+0.42(0.07) 


< +0.28 


|Ni/Fe]„bs 


+0.04(0.05) 


-0.06(0.09) 


-0.02(0.03) 


-0.07(0.07) 


+0.17(0.05) 


< -0.05 


|Ni/Fel„, (EOO) 


+0.03(0.04) 




-0.01(0.02) 


-0.05(0.06) 


+0.09(0.04) 


< -0.04 


|Ni/Fe]„, (Ell) 


+0.04(0.05) 




-0.02(0.03) 


-0.07(0.07) 


+0.16(0.05) 


< -0.05 


|Si/Fe]„bs 


+0.39(0.04) 


+0.43(0.04) 


+0.45(0.02) 


+0.26(0.03) 


+0.40(0.04) 


+0.28(0.08) 


[Si/Fe]c„r (EOO) 


+0.25(0.06) 




+0.18(0.09) 


+0.12(0.10) 


+0.08(0.11) 


+0.16(0.18) 


[Si/Fe]„ (Ell) 


+0.25(0.06) 




+0.18(0.09) 


+0.12(0.10) 


+0.08(0.11) 


+0.16(0.18) 


[Mg/Fe]„bs 


+0.36(0.10) 


+0.43(0.15) 


+0.46(0.05) 


< +0.59 


< +0.29 


< +0.88 


[Mg/Fe]„ (EOO) 


+0.22(0.12) 




+0.19(0.12) 


< +0.45 


< -0.02 


< +0.70 


[Mg/Fe]eor(Ell) 


+0.22(0.12) 




+0.19(0.12) 


< +0.45 


< -0.02 


< +0.70 


|S/Fel„bs 


+0.25(0.08) 


+0.28(0.12) 




+0.27(0.07) 


+0.43(0.09) 


+0.22(0.10) 


|S/Fe]„, (EOO) 


+0.11(0.13) 






+0.13(0.10) 


+0.10(0.18) 


+0.10(0.18) 


|S/Fe],„, (Ell) 


+0.11(0.13) 






+0.13(0.10) 


+0.09(0.18) 


+0.10(0.18) 


|S/Zn|„bs 


+0.02(0.07) 


+0.23(0.15) 




+0.03(0.08) 


+0.00(0.10) 


> -0.06 


|N/Si]„bs 


-0.95(0.07) 


-1.41(0.10) 


> -1.03 


-0.96(0.11) 


-0.98(0.09) 


-1.17(0.10) 


IN/Si]„, (EOO) 


-0.95(0.07) 




> -1.03 


-0.96(0.16) 


-1.00(0.13) 


-1.18(0.10) 


[N/SiU(EIl) 


-0.95(0.07) 




> -1.04 


-0.96(0.16) 


-1.00(0.13) 


-1.18(0.10) 


[N/S]„bs 


-0.81(0.07) 


-1.26(0.14) 




-0.97(0.11) 


-1.01(0.10) 


-1.12(0.10) 



Notes. - The DLA at Zdla - 2.067 toward QB0450-1310 is affected by high ionization effects (see Sect. |2.2K and therefore we do not consider it in this paper. 



- The values in parenthesis coixespond to the 1 tr errors. 

- We provide dust corrections only for the absolute and relative abundances that we will further use in our study. The dust coixections are computed using the lVladilol |20Q2j) prescriptions. The models EOO and El 1 
assume that the intrinsic [Zn/Fe] ratio is equal to +0. 10 dex. 

^ The DLA at Zdla = 2.476 toward QB0841 + 129 does not require any dust conection, because the [Zn/Fe] ratio is almost solar. We obtained for this DLA the rare O abundance measurement which leads to the following 
relative abundances: [0/Fe] = +0.32 ± 0.12, [O/Zn] - +0.26 ± 0.15, and [N/0] = -1.30 ± 0.14. 
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2. 1. Dust depletion corrections 

S everal pieces of evidence show that some dust is present in 
DLAs (iPettiniet al.lll994t iHou et alJlIOOll iProchaska & Wolfd 
12002 '). Three main approaches help to circumvent the problem 
of dust depletion in the studies of DLA abundance patterns. The 
first one, although very limiting, consists in considering only the 
DLAs with low dust depletion levels. [Zn/Fe], the ratio of two 
iron-peak elements with very different degrees of incorporation 
into dust grains, is often used as a dust depletion indicator, with 
Zn being a volatile element and Fe a refractory element. Hence, 
DLAs with [Zn/Fe] < 0.2 are considered as sy stems with negli- 
gible dust depletion levels (' Pettini et al.ll20 00a). This is the case 
in our sample for the DLA at zdla - 2.476 toward Q B0841 + 129 
with [Zn/Fe] - +0.04 + 0.11. In three other DLAs toward 
QB0841-H129, QB1210-H175, and QB2348-1444, the dust de- 
pletion levels are also relatively low: [Zn/Fe] < H-0.28. The 
second approach consists of quantifying the dust depletion ef- 
fects and applying dust corrections to the observed abundances. 
Several authors have considered d ifferent methods to compute 
the dust deple tion corre ctions (e.g .'Vladild'1998US avaglio et alj 
120001; IVladilo 2002a b; IProchas ka & Wolfe 2002). Finally, the 
third approach consists of focusing on non-refractory and/or 
mildly refractory elements, such as N, O, S, Ar, and Zn. In this 
way one assesses the intrinsic abundances in DLAs. 

In our study we considered all the three approaches: the first 
for the DLA atzoLA = 2.476 toward QB0841-H129, and the sec- 
ond and the third for the other DLAs. We performed the dust 
depletio n corrections using the method developed by IVladilol 
(I2002allbh . as in Paper I. This method groups together several 
dust correction models based on different assumptions. We con- 
sider here the models EOO and Ell computed by assuming that 
the intrinsic [Zn/Fe] ratio is equal to -I-O.IO dex (see Vladilo's 
papers for more details). In Table [T] we report the dust-corrected 
values, in addition to the observed ones, for the absolute and rel- 
ative abundances which we will further use in the paper. 



2.2. Ionization corrections 

Several autho r s investigated pho t oioniz at ion in DLAs (e.g . 
ViegasI Il995t iHowk & SembachI 119991: llzotov et al.l l200H 
Vladilo et al.ll200ll) and generally concluded that the ionization 
corrections to observed abundances should be small. We used the 
photoionization diagnostics of Prochaska et al. (2002) to check 
the ionization state of the DLAs in our sample. These diagnostics 
provide a qualitative analysis of the ionization level in a DLA. 
They consist of prescriptions based on the column density ratios 
of adjacent ions of the same element, such as A^(Fe^^)/A^(Fe^), 
and on Ar Ar is very sensitive to ionization, because its pho- 
toionization cross-section is o ne order of magnitude larger than 
for H I at energies hv > 2 Ryd (ISofia & Jenkinsll 19981) . 

These photoionization diagnostics were discussed in Papers I 
and II for the DLAs in our sample. They show that the ionization 
fraction is lower than 10 % for nine out of eleven DLAs and that 
the ionization corrections are negligible for these nine DLAs. 
The two DLAs with the high ionization fractions, the DLA to- 
ward QB2343+125 and the DLA toward QB0450-1310, are 
not considered in the star formation study and therefore the DLA 
toward QB0450-1310 is not included in Table[I] They will be 
discussed in a future paper (Deassuages-Zavadsky et al. in prep). 



3. Chemical evolution models 

Theoretically, we can determine the star formation history and 
the a ge of a galaxy from its abundance patterns (Matteuccil 
l2001b . Indeed, in chemical evolution models the absolute abun- 
dances depend on the model assumptions, whereas the relative 
abundances depend on the nucleosynthesis, the stellar lifetimes, 
and the initial mass function. Relative abundances can thus be 
used as cosmic clocks if they involve two elements formed on 
different timescales. Typical examples are the [a/Fe] and [N/a] 
ratios, since the a-elements (O, Si, S, Mg) are produced on short 
timescales by Type II supernovae (SNe) and the Fe-peak ele- 
ments and nitrogen are produced on long timescales by Type la 
SNe and intermediate-mass stars, respectively. 

In this work, we will make use of different chemical evolu- 
tion models spanning different star formation histories and star 
formatio n efficiencies. Thes e models are the same as the models 
used by Calura et"an (l2003l) and in Paper I. No instantaneous re- 
cycling approximation is adopted, the stellar lifetimes are prop- 
erly taken into account. We describe the models in the next two 
sections. 



3. 1 . Spiral models 

The model for a typical spiral galaxy is calibrated on the chem- 
ical features of the Milky Way. The gal axy is assumed to form 
as a result of two main infall episodes dChiappini et al.|[T997h . 
During the first episode the halo forms and the gas shed by the 
halo rapidly gathers in the center leading to the formation of 
the bulge. During the second episode, a slower infall of exter- 
nal gas gives rise to the disk with the gas accumulating faster 
in the inner than in the outer region, leading to timescales rang- 
ing from ~ 2 Gyr in the inner disk to ~ 7 Gyr in the solar re- 
gion, and up to 13 Gyr in the outer disk ("inside-out" scenario, 
Matteucci & Francois 1989). The efficiency of the star forma- 
tion is characterized by the efficiency parameter, v, defined as 
the inverse of the typical timescale for star formation. It is set 
to y = 1 Gyr"' and becomes zero when the gas surface density 
drops below a certain critical thresho ld. A thre shold density of 
~ 7 M0 pc"^ is adopted in the disk (Chiappini et al. 1997). The 
resulting star formation rate then is proportional to v and the gas 
surface density which depends on the galactocentric radius R and 
time. The initial mass function (IMF) is taken from Scalo ( 198d). 
The galactic wind is assumed to be inefficient ( lMatteuccill200 11) . 

Hence, due to the "inside-out" scenario and the threshold 
density adopted in the spiral model, the star formation history is 
different at different galactocentric radii (see Fig. [T] left panel). 
Indeed, in the inner regions of the disk where the rate of accre- 
tion of matter onto the disk is fast, a high gas surface density well 
above the critical threshold is reached and maintained during a 
long period, and hence the star formation is almost continuous 
through the galaxy lifetime. On the other hand, in the outer re- 
gions of the disk (radius R > 10 kpc) where the rate of accretion 
of matter onto the disk is slow, the SFH proceeds in a gasping 
way, due to the fact that in these regions the gas density is always 
close to the critical threshold. We will label "spiral inner disk" 
the models at /? < 6 kpc, "spiral solar neighborhood" the model 
at /? ~ 8 kpc, and "spiral outer disk" the models at R > 10 kpc. 
These different star formation histories at different galactocen- 
tric radii are associated with different abundance patterns which 
can be compared with observations. The radius R corresponds 
to the position at which the QSO line of sight crosses the DLA 
galaxy disk. 
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Fig. 1. Left panel. Star formation rates per unit area as a function of time of the spiral models at different galactocentric radii. In 
the inner regions of the spiral disk (R < 6 kpc) the star formation proceeds continuously (dotted line), in the solar neighborhood 
(SoN, /? ~ 8 kpc) the star formation also proceeds continuously but at a weaker efficiency (dashed line), and finally in the outer 
regions of the spiral disk {R > 10 kpc) the star formation proceeds in a bursting way (solid line) due to the gas surface density 
threshold adopted for the star formation. Right panel. Star formation rates per unit area as a function of time of the dwarf irregular 
models. Models with two different star formation regimes are illustrated: the dwarf irregular model with a continuous star formation 
(dotted line) and the dwarf irregular model with a bursting star formation (dashed line). The dwarf irregular continuous models are 
characterized by about 10 times weaker star formation rates than the spiral inner disk and solar neighborhood models, and the dwarf 
irregular bursting models are characterized by a single burst of star formation. Here we show a burst lasting 200 Myr and having an 
efficiency of 0.3 Gyr a model which turned out to well represent the star formation history of some DLAs studied in Paper I. 



3.2. Dwarf irregular models 

Dwarf irregular galaxies are assumed to form owing to a contin- 
uous i nfall of pristine gas unti l a mass of ~ 10^ is accumu- 
lated jBradamante et al.lll998 t). The evolution of these galaxies 
is characterized by a bursting star formation history. The param- 
eters which need to be constrained are the number of bursts of 
star formation and for each burst the efficiency, v, (same defini- 
tion as above) and the burst duration, Af . The star formation rate 
is in this case directly proportional to v and the gas density at a 
given time. The star formation in the dwarf irregular models can 
proceed either in short bursts of a duration from 10 to 200 Myr 
separated by long quiescent periods or at a low regime but con- 
tinuously, namely in one or two long episodes of inefficient star 
formation lasting between 3 and 13 Gyr (see Fig.lT] right panel). 
They will be referred as the "dwarf irregular bursting" models 
and the "dwarf irregular continuous" models, respectively. The 
dwarf irregular galaxies are particularly sensitive to outflows re- 
sulting from the ene rgy injection from the supernova explosions 
(lRecchietal.1120011) . 

The IMF is taken from lSalpeterl (Il955l) . The reason for such 
a choice relies mainly on the chemical abundances and on the 
metal content observed in irregulars (for a detailed explana- 
tion, see Calura & Matteucci 2006). Indeed, a Salpeter IMF can 
account for the ab undances observed in local dwarf galaxies 
dRecchi et al.ll2002i ), whereas it leads to an over estimation of the 
abundances in spiral disks (see iRomano et al.] 1 2005) which are 
well accounted for by means of the^ Scalol (119 86) IMF. On the 
other hand, when assuming a Scalo IMF in dwarfs, we underes- 
timate their present-day average metalhcity dCalura & Matteuccil 
l2004l) . 



In summary, three different types of star formation histories are 
modeled by the grid of chemical evolution models which we use 
in our study of DLAs (see Fig. [T]for illustration): 

- the single burst SFH: 

by the dwarf irregular bursting models: dwarf irregular mod- 
els with a single burst of star formation characterized by dif- 
ferent efficiencies v from 0.1 to 5 Gyr ' and different burst 
durations Af from 10 to 200 Myr; 

- the episodic bursting SFH: 

by the spiral outer disk models: spiral models at > 10 kpc 
characterized by a bursting star formation separated by short 
quiescent periods; 

- the continuous SFH: 

by the spiral inner disk models: spiral models at 7? < 6 kpc 
characterized by a continuous star formation over 1 3 Gyr; 
by the spiral solar neighborhood models: spiral models at 
7? ~ 8 kpc characterized by a continuous star formation but 
having a weaker star formation efficiency than the spiral in- 
ner disk models; 

and by the dwarf irregular continuous models: dwarf irreg- 
ular models with a continuous star formation over Af = 13 
Gyr and characterized by low efficiencies v from 0.01 to 0. 1 
Gyr leading to star formation efficiencies about 10 times 
weaker than in the spiral inner disk and solar neighborhood 
models. 

3.3. Stellar yields 

For all the chemical elements studied, the stellar yields con- 
sidered in this work are the same as the ones used in 
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Table 2. Star formation histories, ages, and star formation rates per unit area of the DLA galaxies studied 



QSO 


Zdla 


Possible model 


Star formation characteristics 
v[Gyr-']/Af [Gyr] 


z/ 


Age 
[Gyr] 


SFR per unit area 
[Mo yr-' kpc-2] 


QB0841 + 129 


2.375 


spiral 

dwarf irregular 

dwarf irregular 


outer disk 
bursting: 0.3 / 0.2 

continuous: 0.05 


~ 2.7 
2.6 - 3.0 
-4.0 


0.2 - 0.4 
0.2 - 0.6 
1.0- 1.5 


7.5 X 10-3 
4.0 X 10 3 

1.2 X 10-3 


QB0841 + 129 


2.476 


dwarf irregular 


bursting: 2.0 / 0.02 


-2.5 


0.01 - 0.03 


3.0 X 10 2 


PKS 1157+014 


1.944 


spiral 

dwarf irregular 
dwarf irregular 


outer disk 

bursting: 0.7 /O.l 
continuous: 0.05 


2.0-2.4 
2.4 - 2.9 
4.0 - 5.5 


0.1 - 0.6 
0.7- 1.2 
1.9-2.5 


7.5 X 10-3 

1.6 X 10-2 

7.5 X 10-'' - l.Ox 10-3 






spilal 

dwarf irregular 

dwarf irregular 


illllCl QISK 

bursting: 1.2 / 0.2 

continuous: 0.06 


2.2-2.6 

> 5 


U. i U.J 

0.5 - 1.0 

> 2.4 


O.U A iU i.J A lU 

1.7 x 10-^ 

8.0 X lO-* 


QB2230+02 


1.864 


dwarf irregular 


continuous: 0.1 


>10 


> 3.2 


6.5 X 10" 


QB2348-1444 


2.279 


dwarf irregular 

dwarf irregular 


bursting: 0.1 / 0.15 

continuous: 0.1 


2.4-2.5 

2.6-2.7 


0.1 - 0.3 

0.3 - 0.5 


2.2 X 10-3 

8.7 X lO-* 



Note. The model in boldface corresponds to \he favored model. We assume that the favored model is the model which reproduces 
all the data points within 1 a error (minimum ;^'2) 



Paper I, except for nitrogen. These include the yields of 
iNom oto et'aT' ('1997a') for mass ive stars (M > 10 IVIo), the 
yields of van den Hoeck & Groe nwegenI d 19971) for low- and 
in termediate-mas s stars (0.8 < M/Mq < 8), and the yields 
of INomoto et"an ( f 1997b ) for Type la SNe (model W7). For Zn 
we con sider the specific prescriptions calculated bv lCalura et alj 
(l2003h and for Ni we adopt the prescriptions presented in 
Paper I. 

Our motivation to consider other yields for N comes from 
the lack of reproducing the [N/ff] ratios of the DLAs stud- 
ied in Paper I (see F igs. 15 and 17) and from the work by 
IChiappini et alj (l2003l) . This work shows that by computing the 
same chemical evolution models as those used in this paper 
with the new N yields of Meynet & Maeder (2002) which take 
into account the stellar rotation, the overestimation of the DLA 
[N/a] r atios is solved. The yields with rotation are similar to the 
ones of Ivan den Hoeck & GroenwegenI d 19971) for the low-mass 
stars, but they are substantially smaller for the intermediate-mass 
stars. We henc e adopt the same tables of N yields as the ones of 
IChiappini et a l. (2003) extrapolated from the stellar model grids 
of iMevnet & IVlaeder (200Z) . Meynet & Maeder did not provide 
the yields for other elements yet. The use of a different set of 
yields for N and other elements should, however, not be a prob- 
lem, because the stellar mass range for the production of the 
other elements considered in this work is different. Indeed, O, 
Si, S, and Mg are mainly produced by massive stars, while N is 
mainly produced by intermediate-mass stars. 

Recently. iFranqois et al.l ( |2004|) published a new set of yields 
for several elements studied in this work, in particular Zn, Ni, O, 
and Mg. These yields were derived from a careful comparison 
between the Milky Way chemical evolution model predictions 
and the very accurate UVES measurements of abundance ratios 
observed in Galactic stars with metallicities -4 < [Fe/H] < 0. 
For the elements Ni and Zn, our adopted yields prescriptions (see 
Paper I and Calura et al. 2003, re s pectiv ely) are fully consistent 
with the yields of IFranyois et al.l (|2004|) over the DLA metal- 
licity range (-2.5 < [Fe/H] < -0.5). For O, the new yields 



of iFranqois et al.1 (|2004|) predic t lower [O/Fe] ratios than our 
adopted yields of 'No moto et al.l (Il997ab at a given metallicity. 
This leads to a better agreement with the [O/Fe] measurements 
in DLAs as discussed in Sect. 14.21 The la rgest discrepancy be- 
tween our ad opted yields and the on es of iFranqois et al.l (|2004|) 
concerns Mg. INomoto et al.l (Il997ah yields of Mg are known to 
underestimate the observed [Mg/Fe] ratios in Galactic stars, as 
also shown by Chiappini et al. (1999). We tested the Frangois 
et al. Mg yields, and they usually better reproduce the [Mg/Fe] 
ratios measured in DLAs with the spiral models, but tend to over- 
estimate the [Mg/Fe] ratios with the dwarf irregular models. 



4. DLA star formation histories 

In Paper I we have demonstrated that the star formation history 
and the age of a DLA galaxy can be constrained from a de- 
tailed comparison of the DLA abundance patterns with chemical 
evolution models. The relative abundances when examined as a 
function of a metallicity tracer (e.g. [Fe/H], [Zn/H], or [a/H]) 
allow one to determine the star formation history, and when ex- 
amined as a function of redshift, they determine the age of the 
DLA galaxy defined as the time at which the galaxy started to 
form stars. The star formation rate is also a direct output of the 
model comparison. In the case of the spiral model, the derived 
SFR is the SFR per unit area and it corresponds to the SFR that 
the galaxy has when it is observed. In the case of the dwarf ir- 
regular model, to get the SFR per unit area from the SFR out- 
put, we assume a spherical symmetry and a galactic radius of 
3.5 kpc equivalent to th e mean size of the Magellanic Clouds 
dRussell & Dopitalll992lPl It corresponds to the average SFR in- 
tegrated over the time of the burst of star formation, i.e. over the 
period when the star formation is active. 



3 The radius of dwarf galaxies at high redshift is an unknown quan- 
tity. In the paper, we assume a radius of 3.5 kpc. In the case that the 
radius is half this value, it would produce star fomation rates larger by 
a factor of 4. This difference would, however, have no impact on the 
results described in Sect. [5] 
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The relative abundances that best constrain the chemical 
evolution model are [O/Fe], [Si/Fe], [S/Fe], [S/Zn], [Mg/Fe], 
[Ni/Fe], [N/O], [N/Si], and [N/S], most of them accessible for 
each DLA in our sample (see Table [TJ- For each DLA individu- 
ally, we look for the model that best reproduces all these abun- 
dance ratios as a function of metallicity using the same statistical 
test as in Paper I and defined as follows. We first determine the 
minimal distance between the data point and the curve corre- 
sponding to a given chemical model for each abundance ratio. 
For this purpose, we consider the 1 cr error on the data point 
(i.e. the covariance matrix of the 1 cr measurement error), and 
we derive the minimal distance by computing the distances di 
from the data point to the points defining the considered model 
curve and by determining the di for which the ii,/cr, ratio is min- 
imal. Second, once the minimal distances for all the abundance 
ratios available and for a given model are derived, we compute 
their weighted mean. Finally, the comparison of the weighted 
means obtained for different chemical models determines the 
best chemical evolution model which reproduces the data points 
and thus the DLA galaxy. The limits are not taken into account. 
We then derive the age of the DLA galaxy by determining with 
the same statistical test at which redshift, the so-called redshift of 
formation, the deduced best model reproduces all the abundance 
ratios as a function of redshift. 

We now describe the deduced star formation histories and 
ages of the six DLA galaxies from our new set of data. In Table|2] 
we summarize the results. We provide the best models reproduc- 
ing the DLA abundance patterns and the corresponding derived 
redshifts of formation zy , ages, and star formation rates per unit 
area. 

4.1. QB0841 + 129,Zdi^a=2.375 

The abundance ratios of interest (see above) available in this 
DLA system are [Si/Fe], [S/Fe], [S/Zn], [Mg/Fe], [Ni/Fe], 
[N/Si], and [N/S]. For the refractory elements we consider the 
dust-corrected values given in Table [T] although in this DLA the 
dust corrections are small, since [Zn/Fe] = 4-0.24 + 0.02. 

The spiral model which best reproduces the observed abun- 
dance patterns is the spiral outer disk model (see the thick 
dashed-dotted curve in Fig.|2|. This model fits all the data points 
within 1 o" error, except the [N/Si] ratio which is at 1 .5 cr. We 
also investigated whether a dwarf irregular model can reproduce 
the abundance patterns. We explored both the dwarf irregular 
models with a continuous star formation regime over the Hubble 
time characterized by the efficiencies v = 0.01, 0.03, and 0.05 
Gyr"', and a bursting star formation regime with a burst dura- 
tion At = 0.2 Gyr and characterized by the efficiencies v - 0.1, 
0.3, and 0.5 Gyr '. The bursting model with v = 0.3 Gyr ' re- 
produces all the data points within 1 cr error (thin solid curve in 
Fig.|2]), while the continuous model with v = 0.05 Gyr ' also fits 
well the DLA abundance patterns, except the [S/Zn] ratio which 
is at 1.4 cr (thin dotted curve in Fig.|2l). 

We thus have a degeneracy between these three models 
which all well reproduce the abundance patterns observed in this 
DLA galaxy. However, in terms of derived ages and star forma- 
tion rates, this degeneracy mainly concerns the bursting versus 
continuous star formation regimes. The values given in Table |2l 
indeed, show that the spiral outer disk model and the dwarf irreg- 
ular bursting model, characterized by an episodic bursting and a 
single burst star formation history, respectively (see Sect. [3]), lead 
to very similar results (age and SFR), while those of the dwarf 
irregular continuous model are different. The only way to break 
this degeneracy is to impose that all the abundance patterns have 



to be within 1 cr of the model curve (i.e. have the minimum ;t'^). 
The resulting favored model then is the dwarf irregular burst- 
ing model characterized by v = 0.3 Gyr ' and At = 0.2 Gyr 
However, this selection is very arbitrary. 

From the analysis of the abundance ratios as a function of 
redshift, we obtain a redshift of formation, Zf, of the DLA galaxy 
between 2.6 and 3 for the favored model. This corresponds to an 
age of 0.2 — 0.6 Gyr, the DLA being observed at Zdla — 2.375. 
The star formation rate per unit area is 4 x 10"^ yr"' kpc"^. 

4.2. Q B0841 + 1 29, ZuLA=2.476 

We have at disposal the following abundance ratios in this DLA 
system: [Si/Fe], [S/Fe], [Mg/Fe], [Ni/Fe], [N/Si], [N/S], and the 
rare [O/Fe] and [N/O]. With [Zn/Fe] = 0.04 + 0.11, the system 
is likely free from dust and thus unaffected by dust depletion 
effects. No dust corrections have been computed. 

None of the spiral models succeeds in reproducing the [a/Fe] 
ratios in this DLA galaxy. Indeed, the [O/Fe] ratio is located at 
1.2 cr, [Si/Fe] at 2.4 cr, [S/Fe] at 1.3 cr, and [Mg/Fe] at 1.4 cr 
from the best spiral model curve (see the thick dashed-dotted 
curve in Fig. [3]). Moreover, the [N/O] ratio is also not fitted 
(3.1 cr). Similarly, none of dwarf irregular continuous models 
fits the observed abundance patterns. We investigated dwarf ir- 
regular bursting models with efficiencies between 1.5 to 4 Gyr"^ 
by step of 0.5 Gyr"' and burst durations between 0.01 and 0.05 
Gyr The model with v - 2 Gyr ' and a very short burst duration 
At = 0.02 Gyr best reproduces the observed abundance patterns. 
Indeed, the model fits all the data points within 1 cr error, except 
the [O/Fe] and [N/O] ratios (thin solid curve in Fig.[3]l. 

O is the most reliable diagnostic in the SFH studies, thus the 
fact that the [O/Fe] and [N/O] ratios are not well reproduced by 
the model may suggest either that the O abundance measure- 
ment is underestimated due to a possible O i line saturation (see 
Paper II) or that the adopted O yi elds are not r eliable enough. 
When considering the new yields of|Frangois et al. (2004) which 
predict a lower O abundance at a given metallicity, the problem 
of O is solved for both the dwarf irregular bursting model and 
the best spiral model. 

The degeneracy between the three star formation histories 
investigated is in this case alleviated. Only the dwarf irregular 
bursting model allows to reproduce the abundance patterns of 
this DLA characterized by a high a-enhancement. We would like 
to point out that the adopted Salpeter IMF in the dwarf irregular 
models largely contributes to the successful fit of these DLA data 
points. Indeed, a more top-heavy IMF relative to the Scalo IMF 
is needed to reproduce the high a/Fe ratios. This is why the spiral 
models are ruled out. 

From the abundance ratios versus redshift diagrams, we find 
a redshift of formation Zf - 2.5 for this DLA galaxy, while we 
observe this DLA at zdla - 2.476. This corresponds to a very 
young age of 0.02 + 0.01 Gyr only due to the very short burst du- 
ration, and suggests that the DLA galaxy is likely experiencing 
its first burst of star formation. The measured star formation rate 
per unit area is 3 x 10"^ Mq yr"' kpc"^. Examples of galaxies 
with such a short burst duration exis t among the local metal-poor 
galaxies. In particular. Kunth et alj (Il995h studied the chemical 
properties of IZwl8 by means of a starburst model similar to 
ours, characterized by a burst duration of 0.01 Gyr and an av- 
erage SFR of ~ 0.01 Mq yr"' kpc"^. Transient episodes of star 
formation with such a duration ar e also likely to be ty pical of 
some local blue irregular galaxies (lAnnibah et al.ll2003l) . 
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Fig. 2. Observed and predicted abundance ratios versus metallicity for the DLA at zdla - 2.375 toward QB0841 + 129. The thick 
dashed-dotted curve corresponds to the spiral outer disk model, the thin dotted curve to the dwarf irregular continuous model with 
an efficiency v - 0.05 Gyr and the thin solid curve to the dwarf irregular bursting model characterized by v = 0.3 Gyr ' and a 
burst duration Af = 0.2 Gyr In this and all the following figures, the open squares represent the measured abundan ce ratios, and the 
filled triangles and the filled squares represent the abundance ratios corrected for dust depletion according to the IVladilol ( l2002a[) 
method using the models EGO and Ell, respectively. 



4.3. PKS 7 157+014, Zdi^a=1-944 



Few abundance ratios of interest are available in this DLA sys- 
tem: [Si/Fe], [Si/Zn], [Mg/Fe], [Ni/Fe], and a lower limit to 
[N/Si]. We consider the dust-corrected values given in Table [1] 
since they are relatively important in this system with [Zn/Fe] 
= H-0.37 + 0.07. 



The best spiral model is the spiral outer disk model. It re- 
produces all the data points within 1 cr error and is in agreement 
with the lower limit to [N/Si] (see the thick dashed-dotted curve 
in Fig.|4]i. We considered also dwarf irregular models with both 
the continuous and bursting star formation regimes. For both of 
them we can constrain adequately the star formation parameters 
to reproduce successfully all the abundance patterns within 1 cr 
error: for the continuous model an efficiency v - 0.05 Gyr ' is 
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Fig. 3. Observed and predicted abundance ratios versus metallicity for the DLA at zdla - 2A16 toward QB0841 + 129. The thick 
dashed-dotted curve corresponds to the spiral solar neighborhood model and the thin solid curve to the dwarf irregular bursting 
model characterized by an efficiency v = 2 Gyr ' and a burst duration At = 0.02 Gyr For the definition of symbols, see Fig.|2] 



needed (thin dotted curve in Fig-IH, and for the bursting model 
a burst with an efficiency v = 0.7 Gyr"' and a duration At = 0.1 
Gyr is needed (thin solid curve in Fig. HJl. This leads again to an 
important degeneracy between these three models tracing three 
different star formation histories: an episodic bursting SFH, a 
continuous SFH, and a single burst SFH, respectively. If we 
compare the weighted means of the minimal distances between 
the model curves and the available data points, then the episodic 
bursting SFH (spiral outer disk model) with the minimum is 
the favorite one. However, there is no physical reason to exclude 



the two other star formation histories, since they all agree with 
the abundance patterns. 

From the abundance ratios studied as a function of redshift, 
we derived the redshift of formation and the age of the DLA 
galaxy for the three best models. The obtained age proves to be 
very different from one model to the other Indeed, for the spiral 
outer disk model we find Zf - 2 - 2.4 which corresponds to an 
age of 0.1 - 0.6 Gyr, for the dwarf irregular bursting model we 
have Zf = 2.4 - 2.9 which gives an intermediate age of 0.7 - 1 .2 
Gyr, and for the dwarf irregular continuous model we get Zf = 
4-5.5 which corresponds to a much older age of 1 .9 - 2.5 Gyr. 
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Fig. 4. Observed and predicted abundance ratios versus metallicity for the DLA at zdla - 1 -944 toward PKS 1 157+014. The thick 
dashed-dotted curve corresponds to the spiral outer disk model, the thin dotted curve to the dwarf irregular continuous model having 
an efficiency v = 0.05 Gyr and the thin solid curve to the dwarf irregular bursting model characterized by an efficiency v = 0.7 
Gyr"' and a burst duration Af = 0.1 Gyr. For the definition of symbols, see Fig.|2l 



The star formation rate per unit area is also very disparate from 
one model to the other: 7.5 x 10""' Mq yr"' kpc"^ for the spiral 
model, between 7.5 x 10""* and 1 x lO"-* Mq yr"' kpc"^ for the 
dwarf irregular continuous model, and 1 .6 x 10"^ yr"' kpc"^ 
for the dwarf irregular bursting model. 



4.4. QB1210+175,Zdla=1-892 

We have at disposal the following abundance ratios in this DLA 
system: [Si/Fe], [S/Fe], [S/Zn], [Ni/Fe], [N/Si], [N/S], and an 
upper limit to [Mg/Fe]. We consider the dust-corrected values 
given in Table [T] for the abundance ratios including a refractory 
element, although the dust content in this system is relatively 
low, [Zn/Fe] = +0.24 + 0.07. 

The spiral model which best reproduces the abundance pat- 
terns is the spiral inner disk model. This model fits all the data 
points within 1 cr error, except the [N/S] ratio which is at 2.9 cr 
(see the thick dashed-dotted curve in Fig. |5]). We investigated 
also the dwarf irregular models first to see whether the [N/S] ra- 
tio can better be reproduced and second because the probability 
that the QSO line of sight intercepts the internal parts of a spiral 



galaxy is low. None of the dwarf irregular models with a contin- 
uous star formation regime can reproduce the [N/S] ratio within 
1 cr error. The best of these models is the model with an effi- 
ciency V = 0.06 Gyr" ' , but it fails to fit the [S/Fe] ratio at 1 .4 cr in 
addition to the [N/S] ratio which is at 1.3 cr (thin dotted curve in 
Fig. Is). The favored model is the dwarf irregular bursting model 
with an efficiency v - 1.2 Gyr"' and a duration At = 0.2 Gyr 
(thin solid curve in Fig.|5]l- Indeed, this model reproduces all the 
abundance patterns within 1 cr error and is in agreement with the 
upper limit to [Mg/Fe] as well as the lower limit to [Mg/Fe] ob- 
tained from the Mg ii /12796,2803 doublet which dust-corrected 
value is [Mg/Fe]cor > +0.37 (see Paper II). 

The analysis of the abundance ratios as a function of redshift 
gives for the favored model a redshift of formation z/ = 2.2-2.6 
of this DLA galaxy. This corresponds to an age of 0.5 - 1 Gyr. 
The star formation rate is 1.7 x 10"^ M© yr"' kpc"^. 

4. 5. Q B2230+02, Zqla =1-864 

This DLA galaxy has a very high metallicity, [Zn/H] - -0.67 + 
0.09, higher than 1/5 solar, and relatively important dust de- 
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Fig.5. Observed and predicted abundance ratios versus metallicity for the DLA at zdla - 1.892 toward QB1210+175. The thick 
dashed-dotted curve corresponds to the spiral inner disk model, the thin dotted curve to the dwarf irregular continuous model with 
an efficiency v = 0.06 Gyr and the thin solid curve to the dwarf irregular bursting model characterized by an efficiency v = 1.2 
Gyr ' and a burst duration Af = 0.2 Gyr. For the definition of symbols, see Fig.|2l 



pletion corrections that may bring an uncertainty in the study 
of the star formation history of this object (see Table [T]i. The 
abundance ratios at disposal are [Si/Fe], [S/Fe], [S/Zn], [Ni/Fe], 
[N/Si], [N/S], and an upper limit to [Mg/Fe]. 

Only the spiral inner disk model manages to reproduce the 
observed high metallicity of this DLA galaxy. However, this 
model fails to fit the [N/Si] and [N/S] ratios which are at 2 and 
4.5 cr from the model curve, respectively, and the [Mg/Fe] up- 
per limit (see the thick dashed-dotted curve in Fig. |6]l. We fur- 
ther investigated continuous star formation histories, but with 



weaker star formation efficiencies considering v parameters be- 
tween 0.05 and 0.2 Gyr We found that the dwarf irregular 
continuous model with v = 0.1 Gyr ' best reproduces the abun- 
dance patterns of this DLA galaxy (thin dotted curve in Fig.|6]l. 
It is in agreement with the upper limit to [Mg/Fe] and fits all the 
data points within 1 cr error, except the [N/S] ratio at 2.3 cr. The 
oversolar [Ni/Fe] measurement can hardly be explained by any 
chemical evolution model, since Ni should closely trace Fe. 

We used the abundance ratios versus redshift to estimate the 
redshift of formation and the age of this DLA galaxy for the 
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Fig. 6. Observed and predicted abundance ratios versus metallicity for the DLA at zdla = 1 -864 toward Q B2230+02. The thick 
dashed-dotted curve corresponds to the spiral inner disk model and the thin dotted curve to the dwarf irregular continuous model 
with an efficiency v - 0.1 Gyr For the definition of symbols, see Fig.|2] 



favored model identified. The redshift of formation found is very 
high, Zf > 10, implying a very old age, higher than 3.2 Gyr, for 
this galaxy. This suggests that galaxies were akeady formed at 
z > 10. The star formation rate per unit area of this DLA galaxy 
is of 6.5 X 10-'* Mo yr-' kpc^l 

4.6. Q B2348-1444, Zdla=2.279 

This DLA, on the contrary to the DLA toward Q B2230+02, has 
a low metalHcity 1/180 solar ([Fe/H] = -2.25 ± 0.09). The de- 
rived upper limit to the Zn abundance indicates weak effects due 



to dust depletion, since [Zn/Fe] < +0.28. The abundance ratios 
at disposal for the study of the star formation history of this DLA 
galaxy are [Si/Fe], [S/Fe], [N/Si], [N/S], and the upper limits to 
[Ni/Fe] and [S/Zn]. 

The spiral model which best reproduces the abundance pat- 
terns is the spiral outer disk model. This model fits very well 
the [a/Fe] ratios, but fails to reproduce the [N/Si] and [N/S] ra- 
tios which are at 4.0 cr and 2.3 cr, respectively (see the thick 
dashed-dotted curve in Fig. |7]i. We then considered dwarf ir- 
regular models with both the continuous and bursting star for- 
mation regimes to see whether they better reproduce the data 
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Fig. 7. Observed and predicted abundance ratios versus metallicity for the DLA at zdla - 2.279 toward QB2348-1444. The thick 
dashed-dotted curve corresponds to the spiral outer disk model, the thin dotted curve to the dwarf irregular continuous model with 
an efficiency v = 0.1 Gyr \ and the thin solid curve to the dwarf irregular bursting model characterized by an efficiency v = 0.1 
Gyr"' and a burst duration Af = 0.15 Gyr. For the definition of symbols, see Fig.|2l 



points. The best dwarf irregular continuous model is character- 
ized by an efficiency v = 0.1 Gyr ' (see the thin dotted curve 
in Fig. IT). It reproduces all the abundance patterns within 1 cr, 
except the [Si/Fe] ratio which is at 1 . 1 cr. Similarly, the dwarf 
irregular bursting model, characterized by an efficiency v = 0.1 
Gyr ' and a burst duration Af = 0.15 Gyr, also reproduces all 
the data points within 1 cr, except the [N/Si] ratio which is at 
1.4 cr (see the thin solid curve in Fig. IT). We thus have for this 
DLA galaxy a high degeneracy between these two dwarf irregu- 
lar models tracing two different star formation histories. 



From the abundance ratios studied as a function of redshift, 
we derived the redshift of formation and the age of the DLA 
galaxy for the two possible SFHs. For the dwarf irregular contin- 
uous model we find Zf - 2.6 - 2.7, while for the dwarf irregular 
bursting model we find Zf - 2.4 - 2.5. This corresponds to ages 
of relatively the same order, 0.3 - 0.5 Gyr and 0.1 - 0.3 Gyr, re- 
spectively. The star formation rates per unit area obtained for the 
two models are, on the contrary, relatively different, 8.7 x lO"'* 
and 2.2 x lO^^* yr"' kpc"^, respectively. 



5. Comparison of DLA properties with other 
high-redshift galaxies 

The eleven DLA systems at zdla = 1 .7 - 2.5 studied in Papers I 
and II were selected only on the basis of bright background 
quasars and the existence of HIRES/Keck spectra. They form 
a sample which covers the full range of H i column densities, 
metallicities, and dust depletion levels observed in other DLAs 
at similar redshifts, as illustrated in Fig.[^ However, this sample 
appears to be slightly biased toward high H i column densities, 
since it contains three out of the five higher H i column density 
systems of our control sample, and toward less dusty DLA galax- 
ies, as indirectly selected with the bright QSOs. Nevertheless, 
these effects are not determinant on our results, since they are not 
significant, the full range of DLA values is covered by the DLAs 
in our sample. The derived star formation histories of nine of 



The "control" sample is a compilation of abundance measurements 
from the literature of all DLAs at zdla = 1.7 - 2.5 obtained from high- 
resolution spectra. 



14 



M. Dessauges-Zavadsky et al.: Star formation histories 




20 20,5 21 21.5 

log N(HI) 



-3 -2.5 -2 -1.5 -1 -0.5 
Metallicity [Fe/H] 




0.2 0.4 0.6 0.8 1 
Dust indicator [Zn/Fe] 



Fig. 8. From the left to the right panel. Histograms of the H i column density, metallicity [Fe/H], and dust depletion indicator [Zn/Fe] 
distributions for the whole sample of known DLAs with metallicity measurements in the redshift interval zdla = 1.7 - 2.5 and our 
sample of eleven DLAs studied in Papers I and II (filled area). The comparison of these distributions for the two sets of data shows 
that the DLAs in our sample cover the full range of values for these three properties as observed in DLAs. 



these DLA galaxie^ thus are well representative of the average 
DLA galaxy population. 

Our results show that the abundance patterns observed in 
DLAs are compatible with various star formation histories: 
episodic bursting, single burst, and continuous star formation 
histories with different star formation efficiencies. These are typ- 
ical of inner and outer regions of Milky Way type spiral galaxies 
and of dwarf irregular galaxies with bursting and continuous star 
formation regimes. The degeneracy between these SFHs is un- 
fortunately large for most DLAs. It is only broken for DLAs with 
specific abundance patterns, namely with high a-enhancements, 
high metallicities, and/or low N/a ratios. This is the case of three 
DLAs among the nine in our sample. For the six DLAs with 
degenerated SFHs, we assume the favored SFH is the one for 
which the related model reproduces the abundance patterns with 
the minimum (see Paper I and Table|2l). We can, however, note 
that the main results discussed in the following two sub-sections 
hold despite the degeneracy. 

Let's place the derived DLA properties in the global context 
of high-redshift galaxies and let's compare them with the prop- 
erties of emission-selected galaxies. 

5. 1 . Age distribution 

In Fig.|9](a) we show the age distribution as a function of redshift 
for our nine DLA galaxies (circles), for six emission-selected 
galaxies (stars) which age was derived from the analysis of 
their integrated spectra ( Bruzual 2002|) , and for the Lyman-break 
galaxy (LBG) MS 1512-cB58 (squai'e: [Pettini et alj llOOObl 
[2002). MS 1512-cB58 has age properties representative of the 
LBG population at z ~ 3, estimat ed to be typically between 
40 - 300 Myr (see llwata et ani2005h . 

The DLA galaxies show a large spread in age from a few tens 
of Myr up to a few Gyr The oldest object is the DLA toward 
QB22304-02 with an age larger than 3.2 Gyr. Interestingly, the 
age of this object, observed at zdla = 1 .864, suggests that galax- 
ies were already forming stars at z/ > 10. The same spread in 
age is also observed for the emission-selected galaxies at similar 



^ Two DLAs in our sample have high ionization fractions (see 
Sect. 12. 2t . We exclude them from our star formation history study, be- 
cause the measure of their intrinsic abundances is very uncertain. 



redshifts. We can note that the age distribution of high-redshift 
galaxies is in agreement with the age limit of the Universe de- 
fined by the adopted cosmological model at the epoch of obser- 
vation, but it shows a trend of a steeper age decrease with red- 
shift than the age decrease of the Universe. This could be due to 
a selection bias toward the brightest, strongest star forming and 
hence possibly the youngest emission-selected galaxies at high 
redshift that may not be representative of the galaxy population 
as a whole at z > 1 .5. However, the DLA galaxies selected in ab- 
sorption are independent of such a selection bias, and they also 
show a strong trend toward a very young age of a few hundred 
Myr for galaxies observed at z > 2, much lower than the age of 
the Universe at the epoch of observation. This suggests relatively 
low redshifts of formation (z ~ 3) for most galaxies observed at 
high redshift. 



We investigated different possible coiTelations between the 
age derived for the DLA galaxies in our sample and their Hi 
column density and metallicity properties (see the left panels in 
Fig. [Tol l. No trend is observed as a function of the neutral gas 
content, while we could have expected a decrease of the H i gas 
with the galaxy age increase due to the gas consumption dur- 
ing the star formation. On the other hand, a trend of metallic- 
ity evolution with the age of the DLA galaxies seems to be ob- 
served. However, due to the low number statistics, this correla- 
tion mainly relies on one single measurement made for the DLA 
toward Q B2230-I-02 with the highest age and metallicity. But, 
this measurement is very reliable, since this object is one of the 
rare DLA galaxies where star formation history and age were de- 
termined in a straightforward way; no degeneracy problem was 
encountered (see Sect. 14.5b . Such a coiTelation between age and 
metallicity may simply be explained as the result of the metal 
gas enrichment of a galaxy with time. However, if we consider 
DLA galaxies with the same star formation history, for instance 
only DLAs characterized by a single burst star formation history 
(dwarf irregular bursting models), we can note that the metallic- 
ity is, indeed, on average dependent on the age of the galaxy, but 
it also depends on the star formation efficiency of the burst. This 
is illustrated in the right-bottom panel of Fig.[TO]by the tentative 
correlation between the SFR and metallicity observed for DLAs 
with a single burst SFH. 
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Fig. 9. Panel (a). Age distrib ution as a functi on of redshift for the nine DLA galaxies in our sample (filled circl es), for the emission- 
select ed galaxies studied by iBruzuall (12002 1, stars), and for the Lyman-break galaxy MS 1512-c58 (square; iPettini et al.ll2000R 
I2002h . The solid curve corresponds to the age limit of the Universe defined by the adopted cosmology. Panel (b). Star formation rate 
per unit area distribution as a function of redshift. Our DLA results are shown by filled circles. The dotted-shaded area corresponds 
to the SFRs obtained for the DLAs by Wolfe et al. (2003a) from the Cii* /11335 measurem ents and the dashed- shaded area to 
the SFRs obtained for the emission-selected galaxies from the Gemini Deep Deep Survey (ISavaglio et al ] 12004 . The triangles 
represent the three DLAs for which we have an estimation of their SFR from emission lines (from top to bottom we have the DLAs 
toward QB2206-1958, PKS 0458-02, and HS 1543-H5921, respectively). The square corresponds to the LBG MS 1512-c58 and 
the diamonds to the ORB host galaxies (GRB 000926, ORB 011211, and GRB 021004). Panels (c) and (d). Distributions of the 
star formation histories determined for the DLA galaxies studied as a function of redshift and velocity width of the low-ionization 
absorption lines related to the mass of the galaxy, respectively. Three different SFHs, referred as spiral outer disk (episodic bursting 
SFH), dwarf irregular bursting star formation (single burst SFH), and dwarf irregular continuous star formation (continuous SFH), 
were identified as reproducing the DLA abundance patterns. 



5.2. SXar formation rate distribution 

The star formation rates per unit area derived for the nine DLA 
galaxies in our sample are between -3.2 < logSFR < -1.1 
Mq yr"' kpc"^. These values are i r i very g ood agreement with 
the SFRs obtained bv lWolfeetal] ( l2003ah . -3.0 < log SFR < 
-1.6 Mo yr ' kpc^^, for six DLAs in the redshift interval 



z = n .6. 2.71. IWolfe etal.l ( l2003ah yet used a very different 
technique to determine the DLA SFRs based on the measure 
of the cooling rate of these galaxies assessed by the strength 
of Cii* /11335 absorption. The heating rate is then obtained by 
equating it to the cooling rate, which hence constrains the DLA 
star formation rate, since massive stars forming out of neutral 
gas emit FUV radiation that heats the gas. This shows that the 
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Fig. 10. Investigation for possible correlations between different physical properties - the H i column density, the metallicity, the 
age, and the star formation rate per unit area - derived for the DLA galaxies in our sample. The different symbols correspond to 
the three different star formation histories identified as best reproducing the abundance patterns of a given DLA: spiral outer disk 
(episodic bursting SFH; diamonds), dwarf irregular bursting star formation (single burst SFH; filled circles), and dwarf irregular 
continuous star formation (continuous SFH; crosses). The Kendall rank test was performed to analyze these correlations. We give 
the Kendall rank correlation factor, t, and the probability, P(t), under the null hypothesis of zero correlation (values lower than 5 % 
indicate a significant correlation). 



DLA galaxies seem really to have moderate to low star formation 
rates, on average similar to the SFR per unit area of the Milky 
Way Galaxy estimated to 4 x 10'-^ Mq yr~| kpc"^ dKennicuttl 
119981) . The recent work by I Wild et all (l2007h on Mg ii-selected 
DLAs at, however, lower redshifts 0.4 < zdla < 1 -3 also con- 
firms these results. Indeed, from the detection of the [On] neb- 
ular emission line in 3461 Mg ii absorbers, they derived an aver- 
age SFR per unit area of log SFR » -3.0 M© yr"' kpc"^. 

In Fig. |9](b) we compare the DLA star formation rates per 
unit area (circles) with the SFRs obtained for other high-redshift 
galaxies detected in emission: the star-forming ga laxies from the 
Gem ini Deep Deep Survey (dashed-shaded area; ISavagho et al.l 
120040 . the LBG MS 1512-cB58 (squai-e), and three y-ray burst 



(GRB) host galaxies (diamondsfl The SFR of MS 1512-cB58 
is representative of the average star formati on rate of LBGs 



which have values between 20-370Mb yr~ ' (Pettini et a .119981: 



iKobulnickv & Kooll2000l: llwata et al.ll200l iReddv et al.ll2006^ 



This clearly shows that the DLA galaxies have star formation 
rates several orders of magnitude (10 to 1000 times) lower than 
the SFRs measured in other high-redshift galaxies. Only the 



The GRB 00926 at z = 2.04 l lFvnbo et alj|2002h . the GRB 011211 
at z = 2.14 dFynboetalJ 120031) . and the GRB 021004 at z = 2.34 



dJakobsson et 



ynbo et 
"al.(l2005l) . 
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SFRs of GRB host£| seem to overlap with DLAs with the higher 
SFRs. 

Why are the DLA galaxies so different from other known 
high-redshift galaxies ? There is clearly a strong selection bias 
toward the more luminous and hence toward the stronger 
star forming galaxies in the sample of high-redshift emission- 
selected galaxies. Hence, in a way the emission-selected galax- 
ies are atypical galaxies. The selection of DLA galaxies is, on 
the contrary, largely independent of their luminosity and star 
formation properties. They should thus be representative of the 
whole galaxy population and should account for both the weak 
and strong star forming galaxies. The predominant identifica- 
tion of weak star-forming galaxies with DLAs may hence re- 
flect three indirect effects: (i) high mass and high SFR galaxies 
are rare in the early Universe, therefore a galaxy sample based 
on the covering fraction of the sky (cross-section times number) 
will generally lead to the detection of lower mass and lower SFR 
galaxies; (ii) dust obscuration biases the sample against metal- 
rich DLA galaxies, i.e. against dustier and stronger star forming 
galaxies which could obscure the background quasar (this bias 
is still strongly debated , however; see e.g. Ellison et al. 2001; 
IVladilo & Perouxll2005l) ; and (iii) all galaxies have larger cross- 
section at larger radii, such that DLA sightlines tend to avoid the 
centers of galaxies. This implies that DLA absorptions mostly 
select the external (less active in terms of star formation) regions 
of gala xies as supported by theoretical studies of Mathl in et al.l 
(1200 lb . Although the latter effects certainly contribute, we ex- 
pect the dominant eff'ect is the former: DLAs sample the broad 
distribution of galaxies at high redshift, a distribution dominated 
by low masses and low star formation rates. 

We finally compare the DLA SFRs obtained from absorp- 
tion properties with the few DLAs for which the emission 
counterpart was detected and allowed a direct SFR estima- 
tion (filled triangles in Fig. |9](b)). That conc erns the DLA to- 
ward QB2206-1958 at zdla = 1-919 (e.g. IWeatherlev et all 
I2005h . the DL A SBS 1543-H593 at zpi A = 0.0096 toward 
HS 1543H-5921 (ISchulte-Ladb eck et al.l l2005h. and flie D LA to- 
ward PKS 0458-02 at zdla = 2.039 (M 0llefetan l200l . 
Interestingly, the DLA SFRs inferred from emission show a very 
large dispersion, much larger than the dispersion observed in the 
SFRs derived from absorption properties; with, in particular, the 
DLA galaxy toward QB2206-1958 which reaches a SFR per 
unit area comparable to the median SFR of emission-selected 
galaxies. The question which comes immediately in mind and 
deserves a paper on itself is: is the SFR of the DLA toward 
Q B2206-1958 atypical ? This DLA must in a certain way be ex- 
ceptional, since, despite large observational eff'orts with at least 
20 DLAs at high redshift investigated, emission from the asso- 
ciated galaxies is generally too weak to be detected (see also 
I Wolfe et ani2004 for a discussion). If this DLA is not atyp- 
ical, this would have one major implication: the DLAs would 
hence very likely sample galaxies with a wide range of star for- 
mation efficiencies and may thus cover a much larger part of 
the luminosity function than expected from the SFR distribution 
derived from absor ption only. This is supported by the work of 
IZwaan et al.l (12005.) at z = 0. Based on 21 cm column density 
maps of about 400 galaxies, covering all Hubble types and a 
wide range of luminosities, they showed that the distribution of 
luminosities of the galaxies producing DLAs is nearly flat from 
Mb = -21 to -15. 



' Called also sometimes GRB-DLA galaxies, since they are charac- 
terized by high H i column densities as the DLAs and exhibit damped 
Lycf profiles in their optical afterglow spectra. 



We investigated again different possible correlations be- 
tween this time the star formation rates per unit area obtained 
for the DLA galaxies in our sample and their H i column densi- 
ties and metallicities (s ee the right panels in Fig. [TOli. No clear 
coiTelation is observed. Wolfe et al.l (l2003bl) considered the same 
quantities and also found no correlation between the SFRs and 
the Hi column densities. They argue that such a correlation 
should, however, exist resembling the Kennicutt relation found 
in nearby disk galaxies (Kennicutt 1998). The lack of this corre- 
lation may in our work be due to the fact that with the absorption 
technique we do not sample the global SFR per unit area and the 
global H I column density of the galaxy. The values we get are 
sampled over a transverse dimension corresponding to the linear 
dimension of the quasar, i.e. about 1 pc. A dependence of metal- 
licity on the SFR is also expected. Wolfe et al. found marginal 
evidence for such a coiTelation. We also observe a tentative cor- 
relation when considering only the DLA galaxies with the same 
star formation history, namely those characterized by a single 
burst star formation history (dwarf irregular bursting models). 
This is due to the fact that for these objects we have a more direct 
evidence of star formation feedback leading to higher metallici- 
ties in presence of bursts with higher star formation efficiencies. 

5.3. Star formation history distribution 

Three different star formation histories, referred as spiral outer 
disk (episodic bursting SFH), dwarf iiTegular bursting star for- 
mation (single burst SFH), and dwarf iiTegular continuous star 
formation (continuous SFH), were identified as reproducing the 
DLA abundance patterns (see Sect. [3]l. In Figs. |9](c) and (d), 
we show the distributions of the star formation histories of the 
DLA galaxies studied as a function of redshift and velocity width 
of the low-ionization absorption-lines, respectively. The goal is 
to investigate possible trends of specific star formation histo- 
ries with redshift and the mass of the DLA galaxies. We want, 
however, to underline that both the uncertainty existing for most 
DLAs relative to the degeneracy of their best star formation his- 
tory and the low number of analyzed objects do not allow to draw 
any statistically significant conclusions. The discussion below 
thus is still preliminary. 

The distribution of star formation histories as a function of 
redshift seems to show a trend of finding more galaxies with a 
star formation history typical of dwarf irregulars with a bursting 
star formation regime toward high redshifts, z > 2. This trend 
again supports the first effect mentioned in Sect. 15.21 which im- 
plies that one has a higher probability to encounter younger and 
less evolved galactic structures with lower masses, experienc- 
ing their first bursts of star formation, toward higher redshifts. In 
their study of galactic morphologies in the Hubble Deep Field 
North and South as a function of redshift. [Consehce et alj(l2005l) 
found a clear decline of the number of spirals at z > 1 . It is thus 
not surprising to find only a few DLA galaxies - two out of 
nine - at zdla = 1 .7 - 2.5 having a star formation history typical 
of spiral galaxies. According to the hierarchical theory of galaxy 
formation, these dwarf DLA galaxies at z ~ 2 may represent the 
building blocks constituting the spiral galaxies at z = 0. This 
is also strongly supported by a straightforward mass argument 
which shows that the mass density in DLAs a t z ~ 2 is much 
larger than the mass density in dwarfs today (jProchaska et al.l 
2005). 

We then examined the distribution of star formation histo- 
ries as a function of the velocity width of the low-ionization 
lines, likely related to the mass of DLA galaxies. Indeed, it 
was recently suggested that the observed coiTelation between 
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the metallicity and the velocity width of the low-ionization lines 
in DLAs reveals the existence of a mass-metallicity correlation, 
wit h the velocity width being a mass indicator of these galax- 
ies (IWolfe & ProchaskallT99a: iPeroux et al.ll2003l: iLedoux et alj 
Here, one would expect a trend of finding more galaxies 
with a star formation history typical of spiral outer disks toward 
high masses. No such a trend is observed, the two DLAs with a 
spiral outer disk star formation history do not have particularly 
large velocity widths relative to the six DLAs with a star forma- 
tion history typical of dwarf irregulars with a bursting star forma- 
tion. The DLA galaxy with the largest velocity width has in fact 
a star formation history typical of a dwarf irregular with a con- 
tinuous star formation. The star formation history is in the case 
of this DLA toward Q B2230-I-02 very reliable, since it is one of 
the rare DLAs which star formation history determination is not 
degenerated (see Sect. 14. 51 ). Interestingly, this galaxy is also very 
old with an age longer than 3.2 Gyr As a consequence, the high 
mass of the galaxy may be explained as a requirement to keep 
the star formation on-going during such a long time. 

6. Summary and concluding remarks 

This paper presents the star formation properties of a sample of 
nine DLA galaxies at zdla = 1.7 - 2.5. They were obtained 
thanks to a new technique developed in Paper 1. This technique 
consists of inferring the star formation history and the age of 
DLA systems from a detailed comparison of their intrinsic abun- 
dance patterns (free from dust depletion and ionization effects) 
with chemical evolution models. Comprehensive sets of elemen- 
tal abundances are hence required. Our sample of DLAs fulfills 
this requirement and, we believe, is a sample representative of 
the average DLA galaxy population. Our results are the first step 
toward a possible connection between DLA galaxies and other 
galaxy populations identified in deep imaging surveys with the 
aim of obtaining a global picture of high-redshift objects. We 
obtained the following constraints on the DLA galaxies. 

(1) Are the DLA abundance patterns consistent with a single 
class of star formation history ? The answer is no. The star for- 
mation history was unambiguously constrained for three DLAs, 
and the results clearly show that both a continuous and a sin- 
gle burst star formation history typical of local dwarf irregu- 
lar galaxies are demanded. For the six other DLAs, the star 
formation history is degenerate between these two SFHs, but 
characterized by different star formation burst parameters, and 
the episodic bursting SFH typical of the Milky Way type spi- 
ral outer disks. Among the burst parameters which were in- 
vestigated, only strong star formation bursts with an efliciency 
V >> 5 Gyr ' which predict a far too high a-enhancement at 
low and high metallicities ([Si/Fe] > +0.5 at [Fe/H] < -1) are 
ruled out by the observed DLA abundance patterns. Otherwise, 
the full range of dwarf irregular bursting models with efficien- 
cies V between 0.1 and 4.2 Gyr"' and a relatively short burst 
duration between 0.02 and 0.2 Gyr do satisfy the various DLA 
abundance patterns. Similarly, continuous star formation histo- 
ries characterized by weak star formation efficiencies with v val- 
ues ranging from 0.03 to 0. 1 Gyr ' do agree with the DLA abun- 
dance patterns. However, the spiral inner disk models with a con- 
tinuous star formation history but characterized by a 10 times 
stronger star formation efficiency cannot be ruled out, although 
DLAs intersecting the inner regions of galaxies should be rare. 
Consequently, the abundance patterns show that DLA objects 
do not represent a single class of galaxy population, but sample 
galaxies with a variety of star formation histories. Their common 
characteristic is that all are weak star forming galaxies. 



When adopting the star formation histories related to the 
models which reproduce the abundance patterns with the min- 
imal x^, we do find only two DLA galaxies (each at z < 2) in 
our sample of nine objects with a star formation history typi- 
cal of spiral galaxies. The distribution of the DLA star forma- 
tion histories as a function of redshift shows in fact a trend of 
finding more galaxies with a star formation history typical of 
dwarf irregulars with a bursting star formation toward high red- 
shifts, z > 2. Because DLA galaxies are selected solely on their 
gas cross-section (i.e. independently of their luminosities), this 
trend, if confirmed, indicates, that less evolved galactic struc- 
tures, experiencing their first bursts of star formation and having 
low masses, are more common toward higher redshifts. This is 
further supported by the line profile velocity width, age, and star 
formation rate distributions of the DLA galaxies. 

(2) Are the DLA abundance patterns consistent with a single 
age ? The answer is no again. The DLA galaxies in our sample 
show a large spread in age from a few tens of Myr up to a few 
Gyr. The extremities of the spread, 20 Myr and > 3 .2 Gyr respec- 
tively, correspond to well established ages of the two galaxies 
whose star formation histories are unambiguously constrained. 
DLAs with a-enhanced abundance patterns are associated with 
young objects, since the a-enhancement is a signature of short- 
lived massive stars and thus of a recent burst of star formation. 
Solar a over iron-peak element ratios, in contrast, are associ- 
ated with old objects, due to the time necessary for the release 
of iron-peak elements by Type la SNe and the need of an ineffi- 
cient, continuous star formation. The oldest DLA in our sample 
is observed at zdla = 1 -864 and has an estimated age larger than 
3.2 Gyr It thus nicely indicates that galaxies were already form- 
ing at z/ > 10. 

The general trend which comes out from the DLA age dis- 
tribution as a function of redshift is a much steeper age decrease 
with redshift than the age decrease of the Universe. This is in 
agreement with the age distribution of emission-selected galax- 
ies. A higher proportion of very young galaxies thus seem to 
populate the high redshift Universe at z > 2. This suggests rel- 
atively low redshifts of formation (z ~ 3) for most galaxies ob- 
served at high redshift. 

(3) Are the DLA abundance patterns consistent with a single 
star formation rate ? Although the results indicate a spread of 
two orders of magnitude, this is relatively small in comparison 
with the large range of known star formation rates in galaxies. 
The derived DLA star formation rates per unit area have values 
between -3.2 < logSFR < -1.1 Mp y r '' kpc' ^, in agreement 
with th e SFRs obtained bvl Wolfe et akl (l2003al) and IWild et akl 
(l2007l) using different techniques. In comparison with the SFRs 
of emission-selected galaxies at high redshift, they are 10 to 
1000 times lower Only the SFRs of GRB host galaxies seem 
to overlap with DLAs. 

Why are the DLA galaxies so different from other high- 
redshift objects ? The emission-selected galaxy samples are 
known to be biased toward the more luminous and hence the 
stronger star forming galaxies. But what about DLAs, why 
do they, on the contrary, preferentially select weak star form- 
ing galaxies ? This is likely the result of three indirect effects: 
(i) high mass and high SFR galaxies are rare in the early 
Universe, therefore a galaxy sample based on gas covering- 
fraction of the sky (i.e. the DLA sample), will result in many 
galaxies with lower masses and lower SFRs; (ii) dust obscuration 
may bias the sample against metal-rich and strong star forming 
DLA galaxies; and (iii) the cross-section selection implies that 
DLA sightlines tend to mostly select the external, less active re- 
gions of galaxies. We expect the dominant effect is the former: 
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DLAs sample the broad distribution of galaxies at high redshift, 
and thus this distribution is dominated by young proto-galactic 
structures with low masses and low SFRs. 

At this stage, to improve our knowledge of the star forma- 
tion properties of DLA galaxies, we really need a large sample 
of DLAs with detections of their emission counterparts. Only 
the direct conjunction of absorption and emission properties of 
DLA galaxies will bring a major step in: the comprehension of 
their link with both present-day and other high-redshift galax- 
ies; the determination of their contribution to the star formation 
rate density of the Universe; and the proper probe of the galactic 
chemical evolution. 
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